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^ ' Abstract 

The evolution of (non-strange) antibaryon abundances in the hadronic 
F ^ \ phase of central heavy-ion collisions is studied within a thermal equi- 

ps ' librium framework, based on the well-established picture of subsequent 

■ chemical and thermal freezeout. Due to large annihilation cross sections, 
antiprotons are, a priori, not expected to comply with this scheme. How- 
ever, we show that a significant regeneration of their abundance occurs 

■ upon inclusion of the inverse reaction of multipion fusion, n n ir — > pp 
. (with n n =5-6), necessary to ensure detailed balance. Especially at SPS 

energies, the build-up of large pion-chemical potentials between chemi- 
{Sj ■ cal and thermal freezeout reinforces this mechanism, rendering the p/p 

ratio in reasonable agreement with the observed one (reflecting chemical 
freezeout). Explicit solutions of the pertinent rate equation, which ac- 
count for chemical off-equilibrium effects, corroborate this explanation. 
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: 1 Introduction 

The combined experimental and theoretical study of (ultra-) relativistic heavy- 
ion collisions has lead to impressive progress in our understanding of the mech- 
anisms underlying the formation and evolution of strongly interacting matter. 
A key issue concerns the question whether one is indeed producing systems of 
sufficient macroscopic extent with frequent enough re-interactions among the 
constituents that justify the use of (thermal-) equilibrium techniques. This, 
after all, constitutes an essential prerequisite for a meaningful investigation 
of the phase diagram of QCD in these reactions, including the identification 
of (the approach towards) phase transitions and Quark-Gluon Plasma (QGP) 
formation. A successful interpretation of experimental observables has to be 
based on a comprehensive overall picture of the collision dynamics. For SPS 
energies, this has been accomplished to a remarkable degree, as evidenced by, 
e.g., (i) hadro-chemical equilibration as inferred from produced particle ra- 
tios, (ii) thermal equilibration with collective transverse expansion, as inferred 
from transverse momentum spectra combined with two-particle correlation 



measurements, (iii) a substantial excess of electromagnetic radiation (photons 
and dileptons) ascribed to thermal emission, (iv) the 'anomalous' suppression 
pattern of charmonium states. An important feature of this picture is the 
distinction between chemical and thermal freezeout, implying a significant du- 
ration of an expanding hadronic phase (also essential to explain (ii) and (iii)). 

A solid understanding of the majority of observables then allows to go 
beyond and investigate well-defined deviations. In this talk we will focus on 
antiproton production, which, upon closer inspection, is not easily understood 
within the above picture, but, as we argue below, can be reconciled with it. 

The paper is organized as follows: in Sect. Q we briefly recall the essential 
elements of hadronic equilibration needed to formulate the 'antiproton puzzle'. 
In Sect. [3] we address the problem within a thermal rate equation, first for the 
more transparent equilibrium case, then including off-equilibrium effects via 
its explicit solution. We summarize in Sect. [|. 

2 Equilibration Times and Antiproton Puzzle 

2.1 Two Types of Hadronic Freezeout 

The notion of chemical and thermal freezeout stages in central collisions of 
heavy nuclei can be given a clear physical meaning (and applies equally well to 
hadronic and partonic degrees of freedom). It is based on a hierarchy of time 
scales related to the underlying cross sections. 

Thermal equilibrium in the hadronic fireball can be (locally) maintained as 
long as elastic collisions are able to keep up with its expansion (dilution) rate. 
This criterium can be quantified as 

T 

with Rf a ~ 1/2 (almost independent of collision energy from SIS/BEVALAC 
to SPS) as extracted from hydrodynamic fits to single-particle spectra |l|, 0. 
The thermal relaxation time for a hadron species i is given by 

{Tthd)' 1 = Yl^ihVrel) Ph (2) 
h 

in terms of the density ph of hadrons which serve as scattering centers, and 
the pertinent (thermally averaged) cross section (cr^) (v re i: relative velocity 
between hadron % and h). Standard (resonance dominated) hadronic reac- 
tions, such as 7T7T — > p — > 7T7T, 7iK — > K* — > ttK, ttN — > A — > ttN (or 



even nonresonant iVTV scattering), lead to typical values around (a) ~ 50 mb. 
Even for a rather dilute hadron gas of normal nuclear matter density (ph = 
Po = 0.16 fm -3 ) one finds a thermal equilibration time as short as 2 fm/c. 
A simplified freezeout criterium, which does not explicitly involve the expan- 
sion time scale, is given by the condition that the (thermally averaged) mean 
free path becomes equal to the system size, Xmfp — t fb ( r FB- fireball ra- 
dius). Both criteria, when applied to central heavy-ion collisions at full SPS 
energy (Ei ab = 158 AGeV), indicate thermal freezeout temperatures of about 
T t h = 110 — 120 MeV. This is consistent with combined experimental analysis 
of (light-) hadron pt-spectra and HBT correlations ||, which can be schemat- 
ically summarized by (nonrelativistic) slope parameters T s i ope ^ = T th + rriivf 
with an average transverse flow velocity v t ~ 0.4-0.5 c at thermal freezeout. 

Chemical equilibrium is maintained by inelastic collisions which change 
the particle composition of the system. Corresponding hadronic cross sec- 
tions, e.g., for -Kit «-> KK,pp, are typically much smaller than elastic ones, 
((Ti ne i) ~ 1 mb after thermal averaging. Consequently, chemical equilibration 
times are around T c h e m ~ 100 fm/c, well above possible lifetimes of hadronic 
fireballs created in high-energy nuclear collisions. If anything, the chemical 
composition of the hadronic phase must frozen early in the evolution. At 
full SPS energy essentially all hadron abundances can be accommodated by a 
common chemical freezeout at (T ch , p<£) ~ (170,260) MeV @, §. This result 
is, in fact, difficult to understand in purely hadronic scenarios and has been 
interpreted as 'prehadronic' flavor equilibration, possibly in a QGP. 

A particular consequence of the sequential freezeouts that will be important 
below is the effective conservation of particle numbers for species that are not 
subject to strong decays (e.g., ir, K, r\; after all, the abundances at thermal 
freezeout have to agree with the measured ones). In statistical mechanics 
language this is expressed by finite meson-chemical potentials, which, under 
SPS conditions, reach appreciable values of p 1 ^ ~ 60-80 MeV, p 1 ^ ~ 100- 
130 MeV (similar for p 1 ^), etc. Resonances that are regenerated via strong 
interactions acquire the chemical potentials of their decay products, i.e., p& = 
Pn + Pn, P p = 2p, w , Pk* = Pk + Pw, etc.Q 



1 Strongly decaying resonances with narrow widths (< 50 MeV or so) are in between the 
two cases. However, for the vector mesons lu and <f>, e.g., it turns out that the two limiting 
cases - full chemical equilibrium = 3//^ and /i^ = 2/Ltfc) versus chemical off-equilibrium 
(individual /i w and adjusted to their conserved number) - give numerically very similar 
results. We also note that, if reactions of the type lo~k <-» irir are relevant, the u chemical 
potential would be reduced. 
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Figure 1: Schematic QCD phase dia- 
gram including empirical information 
on chemical freezeout in (ultra-) rela- 
tivistic heavy-ion collisions. 

2.2 Antiproton Systematics 
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Figure 2: Evolution of nucleon- and 
pion-chemical potentials along the 
SPS trajectory of Fig. [TJ, including un- 
certainties in ^ erm ~ 65-80 MeV. 



An enhanced production of antiprotons in central heavy-ion collisions (over the 
extrapolation from pp or pA reactions) has been among suggested signatures 
for QGP formation f?|, |8[] (which is expected to facilitate copious antiquark 
production). The experimental results for p/p ratios at AGS 0, SPS [K| [11 
and RHIC [jT^, [U| are, however, consistent with the chemical freezeout system- 
atics at the respective energies [|, |], H - In Pb(158 AGeV)+Pb, one finds jlO 



p/p = (5.5 ± 1) %, to be compared to p/p = exp(— 2fj, c ^/T c f l ) — 5%. Moreover, 



the p/p ratio does not exhibit a pronounced centrality dependence [fLl]| , which 
is also borne out of hadro-chemical model fits |T5| . 

Despite this apparent agreement it is nevertheless puzzling why the number 
of antiprotons should to be 'frozen in' early in the hadronic phase: large pp 
annihilation cross section of around 50 mb at the relevant thermal energies 
(Ep = m p + §A;T), together with the sizable baryon densities under SPS 
conditions, seem to imply a chemical freezeout for antibaryons at a much lower 
temperature (i.e., later time) than T c h — 170 MeV; e.g., for T = 150 MeV, 
where according to thermal fireball calculations |TB| pb — 0.75 p , one obtains 
for the chemical relaxation time, 



_chem 



{s)v re i) p B 



(3) 



a value of about 3 fm/c, which is significantly smaller than the fireball life- 
time. A naive evaluation of the p/p ratio using thermal freezeout conditions, 
(T th ,/i%) ~ (120,400) MeV pi, yields exp(-2/4/T t7l ) ~ 0.1 %, a factor ~ 50 



below the measured value! This constitutes the 'antiproton puzzle'. 

There have been attempts to resolve this puzzle within sequential scattering 
(transport) simulations [jD], |T5| . The latter are, in principle, suited for a realis- 
tic description of the late (dilute) hadronic stages in heavy-ion collisions. Nev- 
ertheless, extra assumptions, such as a strong in-medium suppression ('shield- 
ing') of the annihilation cross section |T7j or an enhanced production via an 



increased string tension in the early phases |18|], had to be invoked to avoid 
or compensate annihilation losses in the hadron gas. However, a notorious 
problem within these approaches is the treatment of multi-body collisions (3 
or more), which are usually neglected, thus violating detailed balance. For the 
case at hand, these are multi- meson fusion reactions, e.g., n^n — > pp with an 
average of n n ~ 5-6 pions in the incoming channel. In the following we will 
show that the inclusion of these reactions, treated within a thermal rate equa- 
tion, is indeed capable of regenerating a significant number of antiprotons M. 
The notion of detailed balance alone, however, is not enough. A second crucial 
ingredient M is the over-saturation of pion phase-space, encoded in nonzero 
pion-chemical potentials, which enhances the equilibrium abundance of an- 



tiprotons substantially (see also ref. JT9[ for a recent transport approach to the 
problem). 



3 Thermal ^-Production 

3.1 Rate Equation and Chemical Potentials 

In thermal equilibrium, the (net) rate per unit 4-volume for producing antipro- 
tons via n^-pion fusion and its back-reaction, n n ir <-> pp, is given in terms of 
the corresponding invariant matrix element M. nw as (see, e.g., ref. |f20f ) 

Tlf = f d%d%f\d 3 k 5^(K tot ) \M n f {z p z p e"^ - Ce~E~} . (4) 

J i=l 

The thermal factors are in Boltzmann approximation, which allows to factorize 
the fugacity factors z$ = exp(//j/T). The rate equation can then be written in 
the compact form 

n th = dpp^ _PP_ ( _ A 

v dt T? hem \ Z p Zp J ' 

which represents a first order differential equation in the antiproton density pp. 
The relaxation time Tp hem is given in eq. (|3|); the use of the total baryon den- 
sity pb therein (rather than just the proton density p p , which is much smaller) 



(5) 



is based on the assumption that the annihilation cross section of antiprotons 
on other baryons is similar to the one on protons. Eq. (|5|) furthermore pre- 
sumes that p annihilation on excited resonances yields (on average) as many 
additional pions as the resonance itself decays into (e.g., p annihilation on a 
A(1232) gives one pion more than on a nucleon). This leads to a cancellation 
of additional fugacity factors multiplying (z^) and z p in eq. (|5|). The use of 
the total annihilation cross section in Tp hem implies a sum over all possible 
n^-pion final states (ss production is suppressed by the OZI rule). Again, this 
automatically incorporates mesonic resonance states such as p and u to the 
extent that they contribute with a pion-fugacity factor according to their pion 
final states. In a heavy-ion environment this is obviously satisfied for p mesons, 
but to a good approximation also for u (and even 77) mesons, cf. footnote 1. 

Under SPS conditions, where the relative amount of antiprotons in the 
system is small, baryon and pion densities are not significantly altered by the 
reactions under consideration. The rate equation can thus be considered with 
the anitproton density as its only variable. In equilibrium one has IZp = 0, so 
that the antiproton fugacity can be determined as 

zf = z p (z?) . (6) 

Here, as in eq. fl5|), the averaging of z™* is over the initial (or final) pion num- 
bers. Since the pion fugacities enter with large powers, a reasonable estimate 
of Zp requires a reliable description of the final state multiplicities in pp anni- 
hilation, which we will turn to in the following section. 



3.2 Properties of pp Annihilation 

Fig. |3] shows the experimental results for pion multiplicities in pp annihilation 
at rest (see, e.g., the review |H] and references therein). The distribution can 
be reasonably well reproduced by a Gaussian P(n n ) with an average (n^) = 
5.06 and a width Sn^ = 0.9. 

In a thermal heat bath, the annihilation processes occur, on average, at 
somewhat larger CMS energies, ~ 2E p h . Experimental data on the s- 
dependence of (n n ) are displayed in Fig. f|. A linear increase accurately de- 
scribes (n n )(s) in the relevant (thermal) energy range. At the same time, also 
the width of the Gaussian moderately increases. Put together, the averaged 
multi-pion fugacity at given temperature takes the form 

(z?)(T) = J2 p M , (7) 
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Figure 3: Pion-multiplicity distribu- 
tion in pp annihilation at rest, to- 
gether with a Gaussian fit using the 
indicated mean and width. 



10 



A 6 
V 4 



pp -> iyc 




data 
fit 



<n >=2.6+1.3 s 
8n =0.126 <n > s" 



2.5 



3 3.5 
1/2 

s [GeV] 



4.5 



Figure 4: Dependence of average pion- 
multiplicity in pp annihilation on CMS 
energy. The fit and corresponding pa- 
rameterizations are also indicated. 



where for any practical purposes n™ ax = 9. 

For the treatment of off-equilibrium effects (cf. Sect. |3.4[) one also needs 
to account for the energy-dependence of the total annihilation cross section, 
which figures into the equilibration time, eq. @. Its decrease with increasing 
s can be fitted by the empirical formula CTp^ n (s) = (40p ; ~°' 5 + 24 pj ab ' ) mb. 



3.3 p Thermodynamics at SPS Energies 

Inserting the pion-fugacity factor, (z^)(T) from eq. (0) into eq. (|6]), and using 
the temperature-dependent pion- and nucleon-chemical potentials as displayed 
in Fig. [| results in an equilibrium p/p ratio for central Pb(158 AGeV)+Pb 
collisions as represented by the solid lines in Fig. |5|. The enclosed band approx- 
imately reflects the uncertainties in determining the pion-chemical potential, 
associated with the population of higher resonance states (especially those with 
relatively long lifetime) in the late stages of the hadronic evolution. Irrespec- 
tive, the strong enhancement over the equilibrium value without the inclusion 
of finite /i^'s (dashed line in Fig. |5|) is obvious. 

In Fig. |6|, the chemical relaxation time is compared to the fireball lifetime, 
indicating that the antiproton abundance might be driven out of equilibrium 
before thermal freezeout. This will be quantitatively assessed in the following 
section. 
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Figure 5: Equilibrium p/p ratio at Figure 6: Chemical relaxation time for 

SPS energies with zero (dashed line) antiprotons (solid line), eq. @, and 

and finite (solid lines) pion-chemical remaining fireball lifetime until ther- 

potentials (the hatched region corre- mal freezeout (dashed-dotted line) at 



sponds to /C erm =65-80 MeV). 



SPS energies 



3.4 Off-Equilibrium Evolution 

For an explicit solution of the rate equation describing the evolution of the 
antiproton number in an expanding hadron gas we rely on a hierarchy of time 
scales T therm < r^ hem <C r^ em . This implies an equilibrium p density 
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d 3 k 



exp[-(E p - fx p )/T] . 



(8) 



The finite system size in a heavy-ion collision requires the inclusion of volume 
expansion effects [pOp . Using dNp = lZp h VFB(t) dt where Np = PpVfb is the 



antiproton number and Vfb(^) the expanding fireball volume [16|, one obtains 



dpp 
~dt 



-Pp 



+ 



dV, 



FB 



Tp(t) V FB (t) dt 



Tp(t) 



= - Pp (t)L(t) + G(t) (9) 



with gain and loss terms G(t) and L(t), respectively. Eq. (§) is an ordinary 
first order differential equation which has the solution 



Pp {t) = <p{t) ( Pp (t ) + J dt' G(t')Mt') 

to 



exp[- J dt' L(t')] . (10) 



to 



As shown in Fig. [7], the off-equilibrium results confirm the importance of both 
the pion-fugacities as well as the backward reaction of multi-pion fusion in 
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Figure 7: Solid lines: evolution of the p/p ratio according to the solution, 
eq. ( ]T0D , of the rate eq. (0) for an expanding fireball using the free (left panel) 
and a 50% reduced (right panel) o^ n \ long-dashed lines: equilibrium ratio; 
short-dashed lines: solution when switching off the gain term G(t). 

maintaining an antiproton abundance close to the observed value. The effective 
freezeout temperature Tp hem turns out to be close to 130 MeV (slightly higher 
if one assumes a, e.g., 50% smaller in-medium annihilation cross section). 



4 Summary 

We have shown that antiproton production in central heavy-ion collisions at 
SPS energies can be reconciled with the observed value corresponding to stan- 
dard chemical freezeout: annihilation losses in the subsequent hadron gas phase 
are compensated by the back-reaction of multipion fusion, reinforced by the 
build-up of large pion-chemical potentials. We have corroborated our earlier 
findings for the equilibrium case by an explicit solution of the underlying rate 
equation. An interesting test of the regeneration mechanism could be provided 
by the recently suggested balance function technique , which has significant 
sensitivity to the time of (pair) production especially for heavy particles (such 
as nucleons). The off-equilibrium framework is more involved under RHIC 
conditions, as ps is of the same order as p§, so that feedback effects on the 
baryon densities cannot be neglected. This will be addressed in future work, 
as it promises valuable insights into hadro-production at collider energies. 
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